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Abstract 



Apart from omnidirectional, a solid elastic sphere is a natural multi-mode and multi-frequency 
device for the detection of Gravitational Waves (GW). Motion sensing in a spherical GW detector 
thus requires a multiple set of transducers attached to it at suitable locations. If these are resonant 
then they exert a significant back action on the larger sphere and, as a consequence, the joint 
dynamics of the entire system must be properly understood before reliable conclusions can be drawn 
£^ ■ from its readout. In this paper, I present and develop an analytic approach to study such dynamics 

which generalises currently existing ones and clarifies their actual range of validity. In addition, the 
new formalism shows that there actually exist resonator layouts alternative to the highly symmetric 
TIGA, potentially having interesting properties. One of these (I will call it PHC), which only requires 
five resonators per quadrupole mode sensed, and has mode channels, will be described in detail. Also, 
the perturbative nature of the proposed approach makes it very well adapted to systematically assess 
\^ • the consequences of realistic mistunings in the device parameters by robust analytic methods. In 

order to check the real value of the mathematical model, its predictions have been confronted with 
CD . experimental data from the LSU prototype detector TIGA, and agreement between both is found 

to consistently reach a satisfactory precision of four decimal places. 
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! 1 Introduction 

bp. 

The idea of using a solid elastic sphere as a gravitational wave (GW) antenna is almost as old as that 
of using cylindrical bars: as far back as 1971 Forward published a paper Jl6| ] in which he assessed 
some of the potentialities offered by a spherical solid for that purpose. It was however Weber's ongoing 
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philosophy and practice of using bars which eventually prevailed and developed up to the present date, 
with the highly sophisticated and sensitive ultra-cryogenic systems currently in operation — see §§ 



and [11] for reviews and bibliography. With few exceptions [|], spherical detectors fell into oblivion 
for years, but interest in them strongly re-emerged in the early 1990's, and an important number of 
research articles have been published since which address a wide variety of problems in GW spherical 
detector science. At the same time, international collaboration has intensified, and prospects for the 
actual construction of large spherical GW observatories (in the range of ~100 tons) are being currently 
considered in several countries Q even in a variant hollow shape [12]. 

A spherical antenna is obviously omnidirectional but, most important, it is also a natural multi- 
mode device, i.e., when suitably monitored, it can generate information on all the GW amplitudes and 



incidence direction [27], a capability which possesses no other individual GW detector, whether resonant 
or interferometric Furthermore, a spherical antenna could also reveal the eventual existence of 

monopole gravitational radiation, or set thresholds on it ||. 

The theoretical explanation of these facts is to be found in the unique matching between the GW 
amplitude structure and that of the sphere oscillation eigenmodes [24]: a general metric GW generates 



a tidal field of forces in an elastic body which is given in terms of the "electric" components Roioj{t) of 
the Riemann tensor at its centre of mass by the following formula 24]: 

1 There are collaborations in Brazil, Holland, Italy and Spain. 
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f GW (x,t) = £ f^(x) 5 ( /m )(t) (1) 

1=0 and 2 
m=— 



where f (x) are "tidal form factors" , while g^ m ^ (t) are specific linear combinations of the Riemann 
tensor components Roioj(t) which carry all the dynamical information on the GW's monopole (1 = 0) 
and quadrupole (I = 2) amplitudes. It is precisely these amplitudes, g( lm ^ (t), which a GW detector aims 
to measure. 

On the other hand, a free elastic sphere has two families of oscillation eigenmodes, so called toroidal 
and spheroidal modes, and modes within either family group into ascending series of /-pole harmonics, 
each of whose frequencies is (2Z+l)-fold degenerate — see |2j}| for full details. It so happens that 
only monopole and/or quadrupole spheroidal modes can possibly be excited by an incoming metric 
GW and their GW driven amplitudes are directly proportional to the wave amplitudes g^ lm \t) 
of equation (||). It is this very fact which makes of the spherical detector such a natural one for 



GW observations [24|. In addition, a spherical antenna has a significantly higher absorption cross 
section than a cylinder of like fundamental frequency, and also presents good sensitivity at the second 
quadrupole harmonic Q. 

In order to monitor the GW induced deformations of the sphere motion sensors are required. In 
cylindrical bars, current state of the art technology is based upon resonant transducers 0, I7j. A res- 
onant transducer consists in a small (compared to the bar) mechanical device possessing a resonance 
frequency accurately tuned to that of the cylinder. This frequency matching causes back-and-forth res- 
onant energy transfer between the two bodies (bar and resonator), which results in turn in mechanically 
amplified oscillations of the smaller resonator. The philosophy of using resonators for motion sensing 
is directly transplantable to a spherical detector — only a multiple set rather than a single resonator is 
required if its potential capabilities as a multi-mode system are to be exploited to satisfaction. 

A most immediate question in a multiple motion sensor system is: where should the sensors be? 
The answer to this basic question naturally depends on design and purpose criteria. Merkowitz and 
Johnson (M&J) made a very appealing proposal consisting in a set of 6 identical resonators coupling 
to the radial motions of the sphere's surface, and occupying the positions of the centres of the 6 non- 
parallel pentagonal faces of a truncated icosahedron [21, 31]. One of the most remarkable properties 



of such layout is that there exist 5 linear combinations of the resonators' readouts which are directly 
proportional to the 5 quadrupole GW amplitudes g^ 2m \t) of equation (|T|). M&J call these combinations 
mode channels, and they therefore play a fundamental role in GW signal deconvolution in a real, 
noisy system |30|, p5| . In addition, a reduced scale prototype antenna — called TIG A, for Truncated 
Icosahedron Gravitational Antenna — was constructed at Louisiana State University, and its working 



experimentally put to test [29]. The remarkable success of this experiment in almost every detail [32 
|33[ , 34 ] stands as a vivid proof of the practical feasibility of a spherical GW detector H [3(| . 

Despite its success, the theoretical model proposed by M&J to describe the system dynamics is 
based upon a simplifying assumption that the resonators only couple to to the quadrupole vibration 
modes of the sphere [El], |3l| . While this is seen a posteriori of experimental measurements to be a very 



good approximation [29, 33], a deeper physical reason which explains why this happens is missing so 
far. The original motivation for the research I present in this article was to develop a more general 
approach, based on first principles, for the analysis of the resonator problem, very much in the spirit 
of the methodology and results of reference ]^4| ; this, I thought, would not only provide the necessary 
tools for a rigorous analysis of the system dynamics, but also contribute to improve our understanding 
of the physics of the spherical GW detector. 

Pursuing this programme, I succeeded in setting up and solving the equations of motion for the 
coupled system of sphere plus resonators. The most important characteristic of the solution is that it 
is expressible as a perturbative series expansion in ascending powers of the small parameter rj 1 ^ 2 , where 
r] is the ratio between the average resonator's mass and the sphere's mass. The dominant (lowest) 
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order terms in this expansion appear to exactly reproduce Merkowitz and Johnson's equations pl[| , 
whence a quantitative assessment of their degree of accuracy, as well as of the range of validity of their 
underlying hypotheses obtains; if further precision is required then a well defined procedure for going 
to next (higher) order terms is unambiguously prescribed by the system equations. 

Beyond this, though, the simple and elegant algebra which emerges out of the general scheme has 
enabled the exploration of different resonator layouts, alternative to the unique TIG A of M&J. In 



particular, I found one |2^, 26] requiring 5 rather than 6 resonators per quadrupole mode sensed and 



possessing the remarkable property that mode channels can be constructed from the system readouts, 
i.e., five linear combinations of the latter which are directly proportional to the five quadrupole GW 
amplitudes. I called this distribution PHC — see below for full details. 

The intrinsically perturbative nature of the proposed approach makes it also particularly well 
adapted to assess the consequences of small defects in the system structure, such as for example sym- 
metry breaking due to suspension attachments, small resonator mistunings and mislocations, etc. This 
has been successfully applied to account for the reported frequency measurements of the LSU TIGA 
prototype [29], which was diametrically drilled for suspension purposes; in particular, discrepancies be- 
tween measured and calculated values (generally affecting only the fourth decimal place) are precisely 
of the theoretically predicted order of magnitude. 

The method has also been applied to analyse the stability of the spherical detector to several 
mistuned parameters, with the result that it is not very sensitive to small construction errors. This 
conforms again to experimental reports [34], but has the advantage that the argument depends on 
analytic mathematical work rather than on computer simulated procedures — see e.g. [34] or [39]. 

The paper is structured as follows: in section ^, I present the main physical hypotheses of the model, 
and the general equations of motion. In section |3] a Green function approach to solve those equations is 
set up, and in section f| it is used to assess the system response to both monopole and quadrupole GW 



signals. In section 
channels. Section 



5] I describe in detail the PHC layout, including its frequency spectrum and mode 
contains a few brief considerations on the system response to a hammer stroke 
calibration signal, and finally in section [?] I assess how the different parameter mistunings affect the 
detector's behaviour. The paper closes with a summary of conclusions, and three appendices where the 
heavier mathematical details are made precise for the interested reader. 



2 General equations 

With minor improvements, I shall use the notation of references [^] and | |25| |, some of which is now 
briefly recalled. Consider a solid sphere of mass A4, radius R, (uniform) density g, and elastic Lame 
coefficients A and /x, endowed with a set of J resonators of masses M a and resonance frequencies Q a 
(a = l,. . . ,J), respectively. I shall model the latter as point masses attached to one end of a linear 
spring, whose other end is rigidly linked to the sphere at locations x a — see Figure |2[ The system 
degrees of freedom are given by the field of elastic displacements u(x, t) of the sphere plus the discrete 
set of resonator spring deformations z a (t); equations of motion need to be written down for them, of 
course, and this is my next concern in this section. 

I shall assume that the resonators only move radially, and also that Classical Elasticity theory [23] 
is sufficiently accurate for the present purposes^. In these circumstances we thus have 



d 2 



u 



//V 2 u+ (A + /x)V(V-u) +f(x,t) (2.a) 



• dt 2 

z a (t) = -n 2 a [z a (t)-u a (t)]+C xtcrn&l (t) , a=l,...,J (2.b) 



2 We clearly do not expect relativistic motions in extremely small displacements at typical frequencies in the range of 
1 kHz. 



3 




Figure 1: Schematic diagram of the coupling model between a solid sphere and a resonator. The 
notation is that in the text, but sub- indices have been dropped for clarity. The dashed-dotted arc 
line on the left indicates the position of the undeformed sphere's surface, and the solid arc its actual 
position. 



where n a = x a /i? is the outward pointing normal at the the a-th resonator's attachment point, and 

u a {t) = n a -u(x a ,t) , a = l,...,J (3) 

is the radial deformation of the sphere's surface at x a . A dot (') is an abbreviation for time derivative. 
The term in square brackets in ( p.b| ) is thus the spring deformation — q(t) in Figure |2| 

f(x, t) in the rhs of ( |2.a| ) contains the density of all non-internal forces acting on the sphere, which 
is expediently split into a component due the resonators' back action and an external action proper, 
which can be a GW signal, a calibration signal, etc. Then 

f (x, t) = fresonators( X ) t) ~\~ fexternal ( x ; (4) 

Finally, ^ xteinal (t) in the rhs of ([Th]) is the force per unit mass (acceleration) acting on the a-th 
resonator due to external agents. 

Given the hypothesis that the resonators are point masses, the following holds: 

J 

fresonators(x,t) = M a Q? a [z a (t) - U a (t)) (3) (x - X a ) n a (5) 
a=l 

where 5^ is the three dimensional Dirac density function. 

The external forces I shall be considering in this paper will be gravitational wave signals, and also 
a simple calibration signal, a perpendicular hammer stroke. GW driving terms, c.f. equation (|l|), can 
be written 

f GW (x,t) = f( 00 )(x)< 7 ( °)(t) + £ f( 2 ™W 2m )(t) (6) 

m=-2 

for a general metric wave — see [^4| for explicit formulas and technical details. While the spatial 
coefficients f (x) are pure form factors associated to the tidal character of a GW excitation, it is the 
time dependent factors g(' m ) (t) which carry the specific information on the incoming GW. The purpose 
of a GW detector is to determine the latter coefficients on the basis of suitable measurements. 
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If a GW sweeps the observatory then the resonators themselves will also be affected, of course. 
They will be driven, relative to the sphere's centre, by a tidal acceleration which, since they only move 
radially, is given by 

£ w (t) = c 2 R 0i0j {t) x a>i n ajj , a = 1, . . . , J (7) 

where RoiOj(t) are the "electric" components of the GW Riemann tensor at the centre of the sphere. 
These can be easily manipulated to giveQ 

£ W (t) = R £ Y lm (n a )g {lm Ht) , a = l,...,J (8) 

1=0 and 2 

where R is the sphere's radius. 

I shall also be later considering the response of the system to a particular calibration signal, consist- 
ing in a hammer stroke with intensity f"o, delivered perpendicularly to the sphere's surface at point xo: 

fstroke(x, t) = f*o ^ 3 )(x-x )<5(t) (9) 

which is modeled as an impulsive force in both space and time variables. Unlike GW tides, a hammer 
stroke will be applied on the sphere's surface, so it has no direct effect on the resonators. In other 
words, 

e° ke W = 0, a = l,...,J (10) 
The fundamental equations thus finally read: 

d 2 u 

= /^V 2 u + (A + ^)V(V-u) + 

J 

M b n 2 b [z b (t) - u b (t)} 6^(x - x 6 ) n b + f cxtcrnal (x, t) (ll.a) 

6=1 

z a (t) = -n 2 a [z a (t)-u a (t)}+C xter ^\t) , a = l,...,J (ll.b) 

where Internal ( x > t) will be given by either @ or @, as the case may be. Likewise, ^ xternal (i) will be 
given by (||) or (|i~0|), respectively. The remainder of this paper will be concerned with finding solutions 
to the system of coupled differential equations (11), and with their meaning and consequences. 



3 Green function formalism 



In order to solve equations (|ll|) I shall resort to Green function formalism. The essentials of this 
procedure in the context of the present problem can be found in detail in reference 24]; more specific 



technicalities are given in appendix AT| 



By means of such formalism equations ([ll]) become the following integro-differential system: 



u a {t) = uT ternal (t) + J2 Vb f K ab {t-t')[z b {t')-u b {t')]dt' (12.a) 
Za(t) = C Xterna \t) - Q 2 a [za(t) - u a (t)\ , a = 1, . . . , J (12.b) 



3 Yi m (n) are spherical harmonics |15| — see also the multipole expansion of RoiOj{t) in reference |24| 
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where it| xteraal (i) = n a • u external (x a , i) , and u external ( x , t) is the bare (i.e., without attached resonators) 
sphere's response to the external forces f ex ternai( x ) t) in the rhs of Qll.a ). K ab (t) is a kernel matrix 
defined by the following weighted sum of diadic products of wavefunctionsQ: 



K ab (t) = ttlY^ ^n 1 [ n &- u 7v( x &)] [n a -UAr(x a )] smu N t 



N 



Finally, the mass ratios of the resonators to the entire sphere are defined by 



rib 



Mb 
M 



1,...,J 



(13) 



(14) 



and will be small parameters in a real device. 

Before proceeding further, let us briefly pause for a qualitative inspection of equations (0). Equa- 
tion (12.a) shows that the sphere's surface deformations u a (t) are made up of two contributions: one 
due to the action of external agents (GWs or other), contained in u| xternal (i), and another one due to 
coupling to the resonators. The latter is commanded by the small parameters and correlates to 
all of the sphere's spheroidal eigenmodes through the kernel matrix K ab (t). This has consequences for 
GW detectors, for even though GWs may only couple to quadrupole and monopoly spheroidal modes 
of the free sphere [24, Q, attachment of resonators causes, as we see, coupling between these and the 
other modes of the antenna; conversely, the latter back-act on the former, too. As I shall shortly prove, 
such undesirable effects can be minimised by suitably tuning the resonators' frequencies. 



3.1 Laplace transform domain equations 



A solution to equations (12) will now be attempted. Equation (12. a) is an integral equation belonging in 
the general class of Volterra equations jttj, but the usual iterative solution to it by repeated substitution 
of u b (t) into the kernel integral is not viable here due to the dynamical contribution of Zb(t), which is 



in turn governed by the differential equation ( 12. b ). 

A better suited method to solve this integro- differential system is to Laplace-transform it. I denote 
the Laplace transform of a generic function of time f(t) with a caret (") on its symbol, e.g., 



fit) e- st dt 



(15) 



and make the assumption that the system is at rest before an instant of time, t = 0, say, or 



u(x, 0) = u(x, 0) = 2 o (0) = z a (0) = 



(16) 



Equations (12) then adopt the equivalent form 



J 

u a (s) = u^™\ s )-Y / VbKab(s)[z b (s)-u b (s)} (17.a) 

6=1 

s 2 z a (s) = C xtemal (s) - n 2 a [z a (s) - u a (s)] , a = 1, . . . , J (17.b) 

for which use has been made of the convolution theorem for Laplace transforms^. A further simplification 
is accomplished if we consider that we shall in practice be only concerned with the measurable quantities 

4 The capitalised index N will often be used to imply the multiple index {nlm} which characterises the sphere's 
wave-functions. 

° Monopole modes only exist in scalar-tensor theories of gravity, such as e.g. Brans-Dicke |J; General Relativity of 
course does not belong in this category. 

6 This theorem states, it is recalled, that the Laplace transform of the convolution product of two functions is the 
arithmetic product of their respective Laplace transforms. 
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q a {t) = z a (t) - u a (t) , a = l,...,J (18) 

representing the resonators' actual elastic deformations — cf. Figure ^. It is readily seen that these 
verify the following: 



J 



E 

6=1 



S 2 



o~ab + rib 2 Q2 K ab (s) 

S + 5 l n 



„2 ^external I \ 

M = - 7 r^T2 u7 m \s) + ^T^T - « = 1, • • • , J (19) 



Equations (|l9|) constitute a significant simplification of the original problem, as they are a set of 
just J algebraic rather than integral or differential equations. We must solve them for the unknowns 
q a (s), then perform inverse Laplace transforms to revert to q a (t). I come to this next. 

4 System response to a Gravitational Wave 

Our concern now is the actual system response when it is acted upon by an incoming GW. It will be 
calculated by making a number of simplifying assumptions, more precisely: 

) The detector is perfectly spherical. 

) The resonators have identical masses and resonance frequencies. 

) The resonators' frequency is accurately matched to one of the sphere's oscillation eigenfrequencies. 

It will be shown below (section [?]) that a real system can be appropriately treated as one which 
deviates by definite amounts from this idealised construct. Therefore detailed knowledge of the ideal 
system behaviour is essential for all purposes: such is the justification for the above simplifications. 



The wave- functions u n ; m (x) of an elastic sphere can be found in reference [24] in full detail, and I 



shall keep the notation of that paper for them. The Laplace transform of the kernel matrix (|13|) can 



thus be expressed as — see equation (|f8[) in appendix A.l 



Ka b (s) = E -3^T \^(R)\ 2 ^^K-n fc ) = £ -^jrX a f (20) 

nl S +U) nl 47r n l S + ^nl 

where the last term simply defines the quantities Xat>^ ■ Note that the sums here extend over the entire 
spectrum of the solid sphere. 

The assumption that all the resonators are identical simply means that 

m = ... = Vj = rj , fij = . . . = flj = n (21) 

The third hypothesis makes reference to the fundamental idea behind using resonators, which is to 
have them tuned to one of the frequencies of the sphere's spectrum. This is expressed by 

^ = u no i (22) 

where u) no i is a specific and fixed frequency of the spheroidal spectrum. 

In a GW detector it will only make sense to choose Iq = or Iq = 2, as only monopole and quadrupole 
sphere modes couple to the incoming signal; in practice, no will refer to the first, or perhaps second 



harmonic JJj. I shall however keep the generic expression (22) for the time being in order to encompass 
all the possibilities with a unified notation. 

Based on the above hypotheses, equation (^) can be rewritten in the form 
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E 

6=1 



&ab + V E 



2„2 



11/ 



( S 2 + 2 )(s 2 + ^ ) 



(nl) 

2 \ ^afe 



-GW 



s 2 + ft 2 a 



(s) + 



e w (*) 

s 2 + ft 2 



= LOnolo) 



where C? W ( S ) is the Laplace transform of (0), i.e., 



4 GW (*) = ^ E ^K)5 (/m) ( s ) 



(23) 
(24) 



2=0 and 2 
m=—l,...,l 



As mentioned at the end of the previous section, the matrix in the lhs of (23) must now be inverted; 
this will give us an expression for q a (s), whose inverse Laplace transform will take us back to the time 
domain. A simple glance at the equation suffices however to grasp the unsurmountable difficulties of 
accomplishing this analytically. 

Thankfully, though, a perturbative approach is applicable when the masses of the resonators are 
small compared to the mass of the whole sphere, i.e., when the inequality 



7] <C 1 (25) 

holds. I shall henceforth assume that this is the case, as also is with cylindrical bar resonant transducers. 
It is shown in appendix [A.2] that the perturbative series happens in ascending powers of n 1 ^ 2 , rather 
than n itself, and that the lowest order contribution has the form 

?«(*) = r 1 ' 2 E k a m) (s; n) (a) + O(0) , a = 1, . . . , J (26) 

l,m 

where 0(0) stands for terms of order rp or smaller. Here, Aa"^(s; O) is a transfer function matrix which 
relates linearly the system response q a ( s ) to the GW amplitudes g^ lm \s), in the usual sense that q a (t) 
is given by the convolution product of the signal g( lm ' (t) with the time domain expression, 
of A a lm \s;Q). The detector is thus seen to act as a linear filter on the GW signal, whose frequency 
response is characterised by the properties of A a lrn \s;fl). More specifically, the filter has a number of 
characteristic frequencies which correspond to the imaginary parts of the poles of A a lm \s; Q). As also 
shown in appendix |A.2| , these frequencies are the symmetric pairs 



x; l± =tf (l± J^-tl \A nolo (R)\ (a V 1/2 I + 0{n) , a = 1, . . . , J (27) 



47T 

where £ 2 is the a-th eigenvalue of the Legendre matrix 



P lo (n a -n b ) , o,6 = l, J (28) 

associated to the multipole (1$) selected for tuning — see (|^). These frequency pairs correspond to 
beats, typical of resonantly coupled oscillating systems — we shall find them again in section | in a 
particularly illuminating example. 

Equation ( |26| ) neatly displays the amplification coefficient ry -1 / 2 of the resonators' motion ampli- 
tudes, which corresponds to the familiar resonant energy transfer in coupled systems of linear oscilla- 
tors @. 

The specific form of the transfer function matrix Aa (s; O) depends both on the selected mode to 
tune the resonator frequency O and on the resonator distribution geometry. I now come to a discussion 
of these. 
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4.1 Monopole gravitational radiation sensing 

General Relativity, as is well known, forbids monopole GW radiation. More general metric theories, 
e.g. Brans-Dicke ||, do however predict this kind of radiation. It appears that a spherical antenna 
is potentially sensitive to monopole waves, so it can serve the purpose of thresholding, or eventually 
detecting them. It is therefore relevant to consider the system response to scalar waves. 

This clearly requires that the resonator set be tuned to a monopole harmonic of the sphere, i.e., 

O = Lo n0 , (h = 0) (29) 

where n tags the chosen harmonic — most likely the first (n = 1) in a thinkable device. 
Since Pq(z) = 1 (for all z) the eigenvalues of Po( n a are, clearly, 

Cf = J, cf=.--=Cj = o (30) 

for any resonator distribution. The tuned mode frequency thus splits into a single strongly coupled 
pair: 

<4 = V 2 (l± \A n0 (R)\ r, 1 ' 2 j + 0(r/) , n = co n0 (31) 

The A-matrix of equation (^) is seen to be in this case 

AW(s; u> n0 ) = {-If ^ \ [( S 2 + uiy 1 - (s 2 + cu 2 ) 
whence the system response is 



<5«o $m0 (32) 



. / s -1/2 (-!) J 1 



(s'+coiy'-^+^y 



§(«)) (g) + O(0 ) , a = l,...,J (33) 



regardless of resonator positions. The overlap coefficient a n o is given in [p4| p1, and can be calculated by 
means of numerical computer programmes. By way of example, a\o/R = 0.214, and 020/ R = —0.038 for 
the first two harmonics. 

A few interesting facts are displayed by equation (^). First, as already stressed, it is seen that if 
the resonators are tuned to a monopole detector frequency then only monopole wave amplitudes couple 
strongly to the system, even if quadrupole radiation amplitudes are significantly high at the observation 
frequencies uj±. Also, the amplitudes q a (s) are equal for all a, as corresponds to the spherical symmetry 
of monopole sphere's oscillations, and are proportional to J" 1 / 2 , a factor we should indeed expect as 
an indication that GW energy is evenly distributed amongst all the resonators. A single transducer 
suffices to experimentally determine the only monopole GW amplitude g^°°\s), of course, but ( ^3|) 
provides the system response if more than one sensor is mounted on the antenna for whatever reasons. 

4.2 Quadrupole gravitational radiation sensing 

I now consider the more interesting case of quadrupole motion sensing. The choice is now, clearly, 

= ^n2 , (h = 2) (34) 

where n labels the chosen harmonic — most likely the first (n = 1) or the second (n = 2) in a practical 
system. The evaluation of the A-matrix is now considerably more involved l^ql , yet a remarkably 
elegant form is found for it: 



Please note that there is a small notation change: what I call a n o here is a n in 
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k^(s;. n2 ) = (-if J^a n2 E ( E ^ [(s 2 ~ ( 



2 i 2 



6=1 U=^° 



-1 



^ C) ^ C) * 
Cc 



^2m(n b ) (5/2 



(35) 

where Ua is the c-th normalised eigenvector of P2(n a -nft), associated to the non-null eigenvalue Let 
me stress that equation (|35|) explicitly shows that at most 5 pairs of modes, of frequencies cj c ±, couple 
strongly to quadrupole GW amplitudes, no matter how many resonators in excess of 5 are mounted on 
the sphere. The tidal overlap coefficients a n2 can also be calculated (cf. |24|]Q, and give for the first two 
harmonics 



Q12 
R 



0.328 , 



Q22 
R 



0.106 



(36) 



The system response is thus 



Us) = T 1/2 (-i) J \/f ^eIe \ 

X f Yi 
m=-2 

Equation (|37| ) is completely general, i.e., it is valid for any resonator configuration over the sphere's 
surface, and for any number of resonators. It describes precisely how all 5 GW amplitudes g^ 2m \s) 
interact with all 5 strongly coupled system modes; like before, only quadrupole wave amplitudes are seen 
in the detector (to leading order) when 0, = uo n 2 , even if the incoming wave carries significant monopole 
energy at the frequencies uj c ±. 

The degree of generality and algebraic simplicity of ( |37| ) is new in the literature. As we shall now 
see, it makes possible a systematic search for different resonator distributions and their properties. 



,( c )„,( c )* 



n (n b ) g( 2m \s) + O(0) , a = l...,J (37) 



5 The PHC configuration 

Merkowitz and Johnson's TIGA [^] is highly symmetric, and is the minimal set with maximum de- 
generacy, i.e., all the non-null eigenvalues £ a are equal. To accomplish this, however, 6 rather than 5 
resonators are required on the sphere's surface. Since there are just 5 quadrupole GW amplitudes one 



may wonder whether there are alternative layouts with only 5 resonators. Equation (37) is completely 



general, so it can be searched for an answer to this question. In reference [25] a specific proposal was 
made, which I now describe in detail. 

In pursuing a search for 5 resonator sets I found that distributions having a sphere diameter as an 
axis of pentagonal symmetry^ exhibit a rather appealing structure. More specifically, let the resonators 
be located at the spherical positions 

27T 

a = a (all a) , ip a = (a - 1) — , a = 1, . . . , 5 (38) 

5 

The eigenvalues and eigenvectors of i-2( n a -n &) are easily calculated: 



8 Please note that there is a small notation change: what I call a„2 here is b n in p3|. 

9 By this I mean resonators are placed along a parallel of the sphere every 72°. 
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(3 cos 2 a — l) 



ci 2 = e 



15 



sin a cos a 



ci = c 2 



^ sin a 



(m) 



to 



-2, ... ,2 , a = 1, ... ,5 
so the A-matrix is also considerably simple in structure in this case: 



(39.a) 
(39.b) 



where the notation 



4 £ A-l I 
a n2 Sra 2 



2 , 2 



Y 2m (n a )5 l2 , PHC (40) 



u4 ± = O 2 1 ± ^ \A n2 (R)\ Cm r, 1 ' 2 + 0(n) 



m = -2,..., 2 



(41) 



has been used. As seen in these formulas, the five expected pairs of frequencies actually reduce to three, 
so pentagonal distributions keep a certain degree of degeneracy, too. The most important distinguishing 
characteristic of the general pentagonal layout is best displayed by the explicit system response: 



<?a(s) 



-1/2 

V 5 



+ 



+ 



( s 2 + <) _1 -(- 2 + -L) _1 



2Ci 
1 

2(2 



O 9 

s +tu 2+ 



y 1 -{s 2 +ui_) 



Y 20 (n a )g^\s) 

Y 21 (n a )g( u Hs)+Y 2 ^(n a )g( 1 - 1 \ S ) 
Y 22 (n a ) 5 {22) (s) + y 2 -2 K) 5 (2 - 2) (s) 



(42) 



This equation indicates that different wave amplitudes selectively couple to different detector fre- 
quencies. This should be considered a very remarkable fact, for it thence follows that simple inspection 
of the system readout spectrumf*] immediately reveals whether a given wave amplitude g 2m (s) is present 
in the incoming signal or not. 

Pentagonal configurations also admit mode channels, which are easily constructed from (^) thanks 
to the orthonormality property of the eigenvectors ( |39.b ): 



E 



(m)* 



Qa(s) 



r] 1/2 a n2 



2 | 2 
S +UJ m+ 



9 (2m) (s) 



+ O(0) (43) 



These are almost identical to the TIG A mode channels the only difference being that each 
mode channel comes now at a single specific frequency pair co m ±. 



Mode channels are fundamental in signal deconvolution algorithms in noisy systems [3C, 35 1. Pentag- 
onal resonator configurations should thus be considered non-trivial candidates for a real GW detector. 

Based on these facts one may next ask which is a suitable transducer distribution with an axis of 
pentagonal symmetry. Figure ^ shows a plot of the eigenvalues ( |39.a| ) as functions of a, the angular 
distance of the resonator set from the symmetry axis. Several criteria may be adopted to select a 
specific choice in view of this graph. An interesting one can be arrived at by the following argument. If 
for ease of mounting, stability, etc., it is desirable to have the detector milled into a close-to-spherical 



In a noiseless system, of course. 
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polyhedric shaped then polyhedra with axes of pentagonal symmetry must be searched. The number 
of quasi regular convex polyhedra is of course finite — there actually are only 18 of them [|o|, [l8|] — , 
and I found a particularly appealing one in the so called pentagonal hexacontahedron (PHC), displayed 
in Figure ^, left — see also p5fl . This is a 60 face polyhedron, whose faces are the identical irregular 
pentagons of Figure ||, right. The PHC admits an inscribed sphere which is tangent to each face at the 
central point marked in the Figure. It is clearly to this point that a resonator should be attached so as 
to simulate an as perfect as possible spherical distribution. 

The PHC is considerably spherical: the ratio of its volume to that of the inscribed sphere is 1.057, 
which quite favourably compares to the value of 1.153 for the ratio of the circumscribed sphere to the 
TI volume. If the frequency pairs u> m ± are now requested to be as evenly spaced as possible, compatible 
with the PHC face orientations, then the choice a = 67.617° is unambiguously singled out. Hence 



uj ± =ui2 (l± 0.5756 r? 1/2 ) , u 1± = u 12 (l ± 0.8787 ry 1/2 ) , uj 2 ± = u) 12 (l ± 1.0668 r/ 1/2 ) (44) 

for instance for Q = u±2, the first quadrupole harmonic. Figure || shows this frequency spectrum together 
with the multiply degenerate TIG A for comparison. 

The criterion leading to the PHC proposal is of course not unique, and alternatives can be consid- 
ered. For example, if the 5 faces of a regular icosahedron are selected for sensor mounting (a = 63.45°) 
then a four-fold degenerate pair plus a single non-degenerate pair is obtained; if the resonator parallel 
is 50° or 22.6° away from the "north pole" then the three frequencies wo+, + , and u>2+ are equally 
spaced; etc. The number of choices is virtually infinite if the sphere is not milled into a polyhedric 
shape |39|, p§|]. 



Let me finally recall that the complete PHC proposal [25] was made with the idea of building an 
as complete as possible spherical GW antenna, which amounts to making it sensitive at the first two 
quadrupole frequencies and at the first monopole one. This would take advantage of the good sphere 
cross section at the second quadrupole harmonic Q, and would enable measuring (or thresholding) 
eventual monopole GW radiation. Now, the system pattern matrix Aa ^(s;Q) has identical structure 



for all the harmonics of a given I series — see (32) and (p^) — , and so too identical criteria for resonator 



layout design apply to either set of transducers, respectively tuned to u>\2 and u>22- The PHC proposal is 
best described graphically in Figure ||, left: a second set of resonators, tuned to the second quadrupole 
harmonic U22 can be placed in an equivalent position in the "southern hemisphere", and an eleventh 
resonator tuned to the first monopole frequency u>io is added at an arbitrary position. It is not difficult 
to see, by the general methods outlined earlier on in this paper, that cross interaction between these 
three sets of resonators is only second order in i] 1 / 2 , therefore weak. 

A spherical GW detector with such a set of altogether 11 transducers would be a very complete multi- 
mode multi-frequency device with an unprecedented capacity as an individual antenna. Amongst other, 
it would practically enable monitoring of coalescing binary chirp signals by means of a rather robust 



double passage method [10], a prospect which was considered so far possible only with broadband long 
baseline laser interferometers (l^, and is almost unthinkable with currently operating cylindrical 
bars. 



6 A calibration signal: hammer stroke 

This section is a brief digression from the main streamline of the paper. I propose to assess now the 
system response to a particular, but useful, calibration signal: a perpendicular hammer stroke. 



Let us first go back to equation (19) and replace u° xterna (s) in its rhs with that corresponding to a 



hammer stroke, which is easily calculated — cf. appendix [A.l| : 



This is the philosophy suggested and experimentally implemented by Merkowitz and Johnson at LSU. 
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~stroke (s) = _ \A nl (R)\ 2 PKna-no) , a = 1, . . . , J (45) 

nl S +UJ nl 

where no are the spherical coordinates of the hit point on the sphere, and /o ^ no-fo/.M. Clearly, the 
hammer stroke excites all of the sphere's vibration eigenmodes, as it has a completly flat spectrum. 

The coupled system resonances are again those calculated in appendix [A.2| . The same procedures 
described in section || for a GW excitation can now be pursued to obtain 



q a (s) 



-1/2 



(-1) 



j-x 21 + 1 
V 4?r 



fo \A n i(R)\ x 



e(e M^+r 1 



6=1 



2 i 2 



P z (n 6 -n ) + O(0) (46) 



where a = 1,..., J, when the system is tuned to the rji-th spheroidal harmonic, i.e., Vt = u) n i. It is 
immediately seen from here that the system response to this signal when the resonators are tuned to a 
monopole frequency is given by 



? a ( S )=r ? - 1 / 2 (-l) 



sJ-1 



\A n0 (R)\ 



Vt = uj. 



'nO 



(47) 



an expression which holds for all a, and is independent of either the resonator layout or the hit point, 
which in particular prevents any determination of the latter, as obviously expected. The frequencies 
u± are those of (|3l|), and we find here again a global factor J -1 / 2 , as also expected. 

Consider next the situation when quadrupole tuning is implemented, 0, = uj n 2- Only the PHC and 
TIGA configurations will be addressed, as more general cases are not quite as interesting at this point. 



6.1 PHC and TIGA response to a hammer stroke 

Expanding equation ( |46| ) by substitution of the eigenvalues Cm and eigenvectors of the PHC, one 
readily finds that the system response is given by 



Qa{s) = V 



-1/2 f 
JO 



2 

\Aa(R)\ E o (^ + ^ + )" 1 -(^ 2 +^-)' 1 C m 1 Y 2m (n a )Y 2 * m (n ) (48) 



m=-2 



with a = 1, ... ,5, and the mode channels by 



y (m) (s) = v 



-1/2 f 
JO 



A n2 {R)\ 



(s 2 + u 2 m+ ) 1 - ( 



2 i 2 
S +U m _ 



*2m( n 0) 



m 



(49) 



These equations indicate that the system response q a (t) is a superposition of three different beatd^ 
while the mode channels are single beats each, but with differing modulation frequencies. This is 
represented graphically in Figure ||, where we see the result of a numerical simulation of the PHC 
response to a hammer stroke, delivered to the solid at a given location. The readouts q a (t) are somewhat 
complex time series, whose frequency spectrum shows three pairs of peaks — in fact, the lines in the 

ui-)t cos fit, where ui+ and ui- are nearby frequencies, and 

i 10. o 



12 A beat is a modulated oscillation of the form sin | (u>+ 
lu+ +LU- —2Q. The Laplace transform of such function of time is precisely (fi/2) 
higher order terms in the difference lo+ — ui- , which in this case is proportional to r\ 



)-( 



up to 



1/2 
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ideal spectrum of Figure The mode channels on the other hand are pure beats, whose spectra consist 
of the individually separate pairs of the just mentioned peaks. 

The response of the TIGA layout to a hammer stroke has been described in detail by Merkowitz 
and Johnson — see e.g. reference []33|]. The present formalism does of course lead to the results obtained 
by them; in the notation of this paper, we have 



-1/2 



2 i 2 
S +U, 



-1 



2 i 2 
S + LU_ 



-1 



y (m) (s) = -v~ 1/2 fo\A n2 (R)\\ [(, 2 + a;i)" 1 -( S 2 +^) 

for the system response and the mode channels, respectively, where 

T 



P 2 (n a -n ) 
Kjm( n o) , m = -2, 



A n2 {R)\r, 1 l 2 \+0(r ] ) , 



1,...,6 



(50.a) 
,2 (50.b) 

(51) 



are the five-fold degenerate frequency pairs corresponding to the TIGA distribution. Comparison of 
the mode channels shows that they are identical for PHC and TIGA, except that the former come at 
different frequencies depending on the index m. One might perhaps say that the PHC gives rise to 
a sort of "Zeeman splitting" of the TIGA degenerate frequencies, which can be attributed to an axial 
symmetry breaking of that resonator distribution: the PHC mode channels partly split up the otherwise 
degenerate multiplet into its components. 



7 Symmetry defects 

So far we have made the assumption that the sphere is perfectly symmetric, that the resonators are 
identical, that their locations on the sphere's surface are ideally accurate, etc. This is of course un- 
realistic. So I propose to address now how departures from such ideal conditions affect the system 
behaviour. As we shall see, the system is rather robust, in a sense to be made precise shortly, against 
a number of small defects. 

In order to quantitatively assess ideality failures I shall adopt a philosophy which is naturally 
suggested by the results already obtained in an ideal system. It is as follows. 

As seen in previous sections, the solution to the general equations ( |l9| ) must be given as a perturbative 
series expansion in ascending powers of the small quantity r/ 1 / 2 . This is clearly a fact not related to the 
system's symmetries, so it will survive symmetry breakings. It is therefore appropriate to parametrise 
deviations from ideality in terms of suitable powers of n 1 / 2 , in order to address them consistently with 
the order of accuracy of the series solution to the equations of motion. An example will better illustrate 
the situation. 

In a perfectly ideal spherical detector the system frequencies are given by equations (|27|). Now, 
if a small departure from e.g. spherical symmetry is present in the system then we expect that a 
correspondingly small correction to those equations will be required. Which specific correction to the 
formula will actually happen can be qualitatively assessed by a consistency argument: if symmetry 
defects are of order n 1 / 2 then equations (E7|) will be significantly altered in their 77 1 / 2 terms; if on 



the other hand such defects are of order r/ or smaller then any modifications to equations (27) will 
be swallowed into the 0(0) terms, and the more important r/ 1 / 2 terms will remain unaffected by the 
symmetry failure. One can say in this case that the system is robust against that symmetry breaking. 

More generally, this argument can be extended to see that the only system defects standing a chance 
to have any influences on lowest order ideal system behaviour are defects of order 77 1 / 2 relative to an 
ideal configuration. Defects of such order are however not necessarily guaranteed to be significant, and 
a specific analysis is required for each concrete parameter in order to see whether or not the system 



14 



response is robust against the considered parameter deviations. Let us now go into the quantitative 
detail. 

Let P be one of the system parameters, e.g. a sphere frequency, or a resonator mass or location, 
etc. Let -Pideal be the numerical value this parameter has in an ideal detector, and let -P rea i be its value 
in the real case. These two will be assumed to differ by terms of order r, 1 / 2 , i.e., 

Preal = J Pidcal(l+pr/ 1/2 ) (52) 



For a given system, p is readily determined adopting ( |52j ) as the definition of -P rea i, once a suitable 
hypothesis has been made as to which is the value of Pideal- in order for the following procedure to 
make sensible sense it is clearly required that p be of order 1 or, at least, appreciably larger than n 1 ^ 2 . 
Should p thus calculated from (52) happen to be too small, i.e., of order ry 1 / 2 itself or smaller, then the 



system will be considered robust as regards the affected parameter. 
7.1 The suspended sphere 

An earth based observatory obviously requires a suspension mechanism for the large sphere. If a nodal 
point suspension is e.g. selected then a diametral bore has to be drilled across the sphere |29|] . The most 
immediate consequence of this is that spherical symmetry is broken, what in turn results in degeneracy 
lifting of the free spectral frequencies u> n i, which now split up into multiplets uj n i m (m = The 
resonators' frequency f2 cannot therefore be matched to the frequency w no z , but at most to one of 
the 

^nolom's- In this subsection I keep the hypothesis — to be relaxed later, see below — that all the 
resonators are identical, and assume that Q falls within the span of the multiplet of the t(; no z om 's. Then 

Un l m = ^ 2 (! +PmV 1/2 ) , m = -l , . . . ,l (53) 



The coupled frequencies, i.e., the roots of equation (p6|), will now be searched. The kernel matrix 
K a b(s) is however no longer given by (p0[), due the removed degeneracy of uj n i, and we must stick to its 
general expression (|75"|), or 



«M = E \ A m(R)\ 2 ^ iWna) - E iar^ x ^ m) (54) 

nlm S +UJ nlm 47r n / m S + ^nlm 



Following the steps of appendix A.l we now seek the roots of the equation 



det 



+ Tj 



n 2 s 2 



(nolom) 
Xab +V 



E 



n 2 s 2 

; s 2 + 2)( s 2 + 2 



(nlm) 
Xab 







m=-l ( s2 + ^ 2 )( s2+U; noZom) „/---„„/„.,„ ^ ^" ' 1 " ^ ^ "hn > 

(55) 

Since O relates to uj no i m through equation ( |53| ) we see that the roots of ( |55| ) fall again into either of 
the two categories (|8^) (see Appendix A. 2 ), i.e., roots close to ±zf2 and roots close to ±iuj n i m (nl ^ uqIq). 
I shall exclusively concentrate on the former now. Direct substitution of the series (fOg ) into (H) yields 
the following equation for the coefficient xi : 



det 



Sab 



(n l m) 
Xab 



XI 



Pr, 







(56) 



This is a variation of ([K]), to which it reduces when p m = 0, i.e., when there is full degeneracy. 
The solutions to ( |56| ) no longer come in symmetric pairs, like (|27|). Rather, there arc 21Q+1+J of 
them, with a maximum number of 2(2/o+l) non-identically zero roots if J > 2/p+lp^|. For example, if 

13 This is a mathematical fact, whose proof is relatively cumbersome, and will be omitted here; let me just mention that 
it has its origin in the linear dependence of more than 2?o+l spherical harmonics of order Zq. 
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we choose to select the resonators' frequency close to a quadrupole multiplet (Zq = 2) then (^) has at 
most 5+ J non-null roots, with a maximum ten no matter how many resonators in excess of 5 we attach 
to the sphere. Modes associated to null roots of ( |56| ) can be seen to be weakly coupled, just like in a 
free sphere, i.e., their amplitudes are smaller than those of the strongly coupled ones by factors of order 

In order to assess the reliability of this method I have applied it to see what are its predictions for 
a real system. To this end, data taken with the TIGA prototype at LSl/[^ were used to confront with. 
The TIGA was drilled and suspended from the centre, so its first quadrupole frequency split up into a 
multiplet of five frequencies. Their reportedly measured values are 



U120 = 3249 Hz , wiai = 3238 Hz , uj 12 _i = 3236 Hz , w 122 = 32 24 Hz , uj 12 _ 2 = 3223 Hz , (57) 
All 6 resonators were equal, and had the following characteristic frequency and mass, respectively: 

n = 3241 Hz , n = (58) 

' 1762.45 v ; 

Substituting these values into (|53|) it is seen that 

Po = 0.2075 , pi = -0.0777 , p_i = -0.1036 , p 2 = -0.4393 , p_ 2 = -0.4650 (59) 



Equation (56) can now be readily solved, once the resonator positions are fed into the matrices 
Xab • Such positions correspond to the pentagonal faces of a truncated icosahedron. Merkowitz [29] 



gives a complete account of all the measured system frequencies as resonators are progressively at- 
tached to the selected faces, beginning with one and ending with six. Figure ^ graphically displays the 
experimentally reported frequencies along with those calculated theoretically by solving equation (|5q). 
In Table 1 I give the numerical values. As can be seen, coincidence between theoretical predictions and 
experimental data is remarkable: the worst error is 0.2%, while for the most part it is below 0.1%. This 
is a few parts in 10 4 , which is precisely the magnitude of rj, as specified in equation (|5S|). 

Therefore discrepancies between theoretical predictions and experimental data are exactly as expected, 
i.e., of order rj. In addition, it is also reported in reference |33| that the 11-th, weakly coupled mode of 
the TIGA (enclosed in square brackets in Table 1) has a practically zero amplitude, again in excellent 
agreement with the general theoretical predictions about modes beyond the tenth — see paragraph after 



equation (56) 



This is a remarkable result which encouraged a better fit by estimates of next order corrections, i.e., 



Xi of ( 89. a ). As it turned out, however, matching between theory and experiment does not consistently 



improve in the next step. This is not really that surprising, though, as M&J explicitly state [33] that 



control of the general experimental conditions in which data were obtained had a certain degree of 
tolerance, and they actually show satisfaction that ~1% coincidence between theory and measurement 
is comfortably accomplished. But 1% is two orders of magnitude larger than r\ — cf. equation ( |58"1) — , 
so failure to refine our frequency estimates to order r\ is again fully consistent with the accuracy of 
available real data. 

A word on a technical issue is in order. Merkowitz and Johnson's equations for the TIGA |2l] , |3l| 
are identical to the equations in this paper to lowest order in r/. Remarkably, though, their reported 
theoretical estimates of the system frequencies are not quite as accurate as those in Table 1 j33|. The 
reason is probably this: in M&J's model these frequencies appear within an algebraic system of 5+J 
linear equations with as many unknowns which has to be solved; here instead the algebraic system has 



only J equations and unknowns, actually equations (IS). This is a very appreciable difference for the 
range of values of J under consideration. While the roots for the frequencies mathematically coincide in 
both approaches, in actual practice they are estimated, generally by means of computer programmes. 



These data are contained in reference |E9| ; I want to thank Stephen Merkowitz for kindly handing them to me. 
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Table 1: Numerical values of measured and theoretically predicted frequencies (in Hz) for the TIG A 
prototype with varying number of resonators. Relative errors are also shown as parts in 10 4 . The 
calculated values of the tuning and free multiplet frequencies are taken by definition equal to the 
measured ones, and quoted in brackets. In square brackets the frequency of the weakly coupled sixth 
mode in the full, 6 resonator TIGA layout. These data are plotted in Figure |(| 



Descr. 


Meas. 


Calc. 


Difference 
(parts in 10 4 ) 


Descr. 


Meas. 


Calc. 


Difference 
(parts in 10 4 ) 


Tuning 


3241 


(3241) 


(0) 


4 reson. 


3159 


3155 


-12 


No reson. 


3223 


(3223) 


(0) 




3160 


3156 


-11 




3224 


(3224) 


(o) 




3168 


3165 


-12 




3236 


(3236) 


(0) 




3199 


3198 


-5 




3238 


(3238) 


(0) 




3236 


3236 







3249 


(3249) 


(0) 




3285 


3286 


3 


1 reson. 


3167 


3164 


-8 




3310 


3310 







3223 


3223 







3311 


3311 







3236 


3235 


-2 




3319 


3319 







3238 


3237 


-2 


5 reson. 


3152 


3154 


8 




3245 


3245 







3160 


3156 


-14 




3305 


3307 


6 




3163 


3162 


-3 


2 reson. 


3160 


3156 


-13 




3169 


3167 


-8 




3177 


3175 


-7 




3209 


3208 


-2 




3233 


3233 







3268 


3271 


10 




3236 


3236 







3304 


3310 


17 




3240 


3240 







3310 


3311 


3 




3302 


3303 


3 




3313 


3316 


10 




3311 


3311 







3319 


3321 


6 


3 reson. 


3160 


3155 


-15 


6 reson. 


3151 


3154 


11 




3160 


3156 


-13 




3156 


3155 


-3 




3191 


3190 


-2 




3162 


3162 







3236 


3235 


-2 




3167 


3162 


-14 




3236 


3236 







3170 


3168 


-7 




3297 


3299 


8 




[3239] 


[3241] 


[6] 




3310 


3311 


2 




3302 


3309 


23 




3311 


3311 







3308 


3310 


6 












3312 


3316 


12 












3316 


3317 


2 












3319 


3322 


10 
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It is here that problems most likely arise, for the numerical reliability of an algorithm to solve matrix 
equations normally decreases as the rank of the matrix increases. 

7.2 Other mismatched parameters 

We now assess the system sensitivity to small mismatches in resonators' masses, locations and frequen- 
cies. 

7.2.1 Resonator mass mismatches 

If the masses are slightly non-equal then one can write 

M a = V M(l + fiaV 1/2 ) , a = l,...,J (60) 

where n can be defined e.g. as the ratio of the average resonator mass to the sphere's mass. It is 
immediately obvious from equation (|60|) that mass non-uniformities of the resonators only affect the 
main equations in second order, since resonator mass non-uniformities result, as we see, in corrections 
of order ry 1//2 to r/ 1 / 2 itself, which is the very parameter of the perturbative expansions. The system is 
thus clearly robust to mismatches in the resonator masses of the type (|60|). 



7.2.2 Errors in resonator locations 

The same happens if the locations of the resonators have tolerances relative to a preselected distribution. 
For let n a be a set of resonator locations, for example the TIGA or the PHC positions, and let n' a be 
the real ones, close to the former: 

n' a = n a + v a rj 1/2 , a = 1, . . . , J (61) 



The values n a determine the eigenvalues £ a in equation (27), and they also appear as arguments 
to the spherical harmonics in the system response functions of sections [|-|6]. It follows from ( |6l|) by 
continuity arguments that 



Yi m (n' a ) 

c 



Y lm (n a ) + Oiv 1 / 2 ) 



(62.a) 
(62.b) 



Inspection of the equations of sections shows that both Q a and Yj m (n a ) always appear within 



lowest order terms, and hence that corrections to them of the type (62) will affect those terms in second 
order again. We thus conclude that the system is also robust to small misalignments of the resonators 
relative to pre-established positions. 

7.2.3 Resonator frequency mistunings 

The resonator frequencies may also differ amongst them, so let 



Q. a = Q (1 + p a rj 



1/21 



1,...,J 



(63) 



To assess the consequences of this, however, we must go back to equation (86) and see what the 
coefficients in its series solutions of the type ( |89.a| ) are. The procedure is very similar to that of 
section 7.1, and will not be repeated here; the lowest order coefficient xi is seen to satisfy the algebraic 

2 

equation 



det 



c 1 Xac Ocb 

lab 2^ 

XI ±^ Xl ~ Pc 

9 C— U o 



(64) 
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which reduces to (|90| ) when all the p's vanish, as expected. This appears to potentially have significant 
effects on our results to lowest order in ?? 1//2 , but a more careful consideration of the facts shows that 
it is probably unrealistic to think of such large tolerances in resonator manufacturing as implied by 
equation (^3|) in the first place. In the TIGA experiment, for example J29|, an error of order r/ 1 / 2 would 
amount to around 50 Hz of mistuning between resonators, an absurd figure by all means. In a full scale 
sphere (~40 tons, ~3 metres in diameter, ~800 Hz fundamental quadrupole frequency, r/~10 -5 ) the 
same error would amount to between 5 Hz and 10 Hz in resonator mistunings for the lowest frequency. 
This is probably excessive for a capacitive transducer, but may be realistic for an inductive one. With 
this exception, it is thus more appropriate to consider that resonator mistunings are at least of order 
rj. If this is the case, though, we see once more that the system is quite insensitive to such mistunings. 

Summing up the results of this section, one can say that the resonator system dynamics is quite 
robust to small (of order r/ 1//2 ) changes in its various parameters. The important exception is of course 
the effect of suspension drilling, which do result in significant changes relative to the ideally perfect 
device, but which can be relatively easily calculated. The theoretical picture is fully supported by 
experiment, as robustness in the parameters here considered has been reported in the real device |33| . 



8 Conclusions 

A spherical GW antenna is a natural multi-mode device with very rich potential capabilities to detect 
GWs on earth. But such detector is not just a bare sphere, it requires a set of motion sensors to be 
practically useful. It appears that transducers of the resonant type are the best suited ones for an 
efficient performance of the detector. Resonators however significantly interact with the sphere, and 
they affect in particular its frequency spectrum and vibration modes in a specific fashion, which must 
be properly understood before reliable conclusions can be drawn from the system readout. 

The main objective of this paper has been the construction and development of an elaborate theoret- 
ical model to describe the joint dynamics of a solid elastic sphere and a set of radial motion resonators 
attached to its surface at arbitrary locations, with the purpose to make predictions of the system 
characteristics and response, in principle with arbitrary mathematical precision. 

The solutions to the equations of motion have been shown to be expressible as an ascending series 
in powers of the small "coupling constant" r], the ratio of the average resonator mass to the mass of the 
larger sphere. The lowest order approximation corresponds to terms of order r/ 1 / 2 and, to this order, 
previous results fffi , B9 , |2^| are recovered. This, I hope, should contribute to clarify the nature of the 
approximations inherent in earlier approaches, and to better understand the physical reason for their 
remarkable accuracy [j33j ]. 

In addition, the methods of this paper have permitted us to discover that there can be in fact trans- 
ducer layouts alternative to the highly symmetric TIGA, and having potentially interesting practical 
properties. An example is the PHC distribution, which is based on a pentagonally symmetric set of 5 
rather than 6 resonators per quadruopole mode sensed. This transducer distribution has the property 
that mode channels can be constructed from the resonators' readouts, much in the same way as in the 



TIGA [31|. In the PHC however a new and distinctive characteristic is present: different wave ampli- 
tudes selectively couple to different detector modes having different frequencies, so that the antenna's 
mode channels come at different rather than equal frequencies. The PHC philosophy can be extended 
to make a multi-frequency system by using resonators tuned to the first two quadrupole harmonics of 
the sphere and to the first monopole, an altogether 11 transducer set. 

The assessment of symmetry failure effects, as well as other parameter departures form ideality, 
has also been subjected to analysis. The general scheme is again seen to be very well suited for the 
purpose, as the theory transparently shows that the system is robust against relative disturbances of 
order n or smaller in any system parameters, also providing a systematic procedure to assess larger 
tolerances — up to order r/ 1 / 2 . The system is shown to still be robust to tolerances of this order in 
some of its parameters, whilst it is not to others. Included in the latter group is the effect of spherical 
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symmetry breaking due to system suspension in the laboratory, which causes degeneracy lifting of the 
sphere's eigenfrequencies, which split up into multiplets. A strong point is that, by use of mostly 
analytic algorithms, it has been possible to accurately reproduce the reportedly measured frequencies 
of the LSU prototype antenna [E9] with the predicted precision of four decimal places. The also 



reported robustness of the system to resonator mislocations [ 33 [ is too in satisfactory agreement with 
the theoretical predictions. 

The perturbative approach here adopted is naturally open to refined analysis of the system re- 
sponse in higher orders in rj. For example, one can systematically address the weaker coupling of 
non-quadrupole modes, etc. It appears however that such refinements will be largely masked by noise 



in a real system, as shown by Merkowitz and Johnson [34], and this must therefore be considered first. 
So the next step is to include noise in the model and see its effect. Stevenson [39] has already made some 
progress in this direction, and partly assessed the characteristics of TIGA and PHC, but more needs 
to be done since not too high signal-to-noise ratios should realistically be be considered in an actual 
GW detector. In particular, mode channels are at the basis of noise correlations and dependencies, as 



well as the errors in GW parameter estimation [35]. I do expect the analytic tools developed in this 



article to provide a powerful framework to address the fundamental problems of noise in a spherical 
GW antenna which, to my knowledge, have not yet received the detailed attention they require. 
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Appendices 

A.l Green functions for the multiple resonator system 

The density of forces in he rhs of equation ( |2.a| ) happens to be of the separable type 

f(x,t)=£f (a) (x)0 (a) (t) (65) 



where a is a suitable label. It is recalled from reference [24] that, in such circumstances, a formal 
solution can be written down for equation (|2.a| ) in terms of a Green function integral, whereby the 
following orthogonal series expansion obtains: 



u(x, *) = ££ ^ f^ u N { X ) gP (t) (66) 

a N 

where 



# = ^7 / u^W-f (a) (x)A (67.a) 



9P(t) = /V Q) (Osin^-0^' (67.b) 



o 



Here, ujn and uat(x) are the eigenfrequencies and associated normalised wave- functions of the free 
sphere. Also, N is an abbreviation for a multiple index {nlm}. The generic index a is a label for the 
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different pieces of interaction happening in the system. I quote the result of the calculations of the 
terms needed in this paper: 



/relnators.TV = fi a [ n a"U^(x a )] , a = 1, . . . , J (68.a) 

/gw,jv = a m5wbmm> , N = {nlm} , l' = 0, 2 , m' = -/',..., I' (68.b) 

/stroke.TV = M^ 1 f -U^(x ) (68.c) 

where the coefficients a n i in ( |68.b ) are overlapping integrals of the type ( 67.a ), and 



Resonators, jv(*) = / WO - «o(*)] sina;jv(i - dt' , a = 1, . . . , J (69.a) 
<?GWiv(*) = f 9 {lm) {t') S \nu: N {t-t')dt' (69.b) 



3stroke,iv(i) = sinwArt (69.c) 
If this is replaced into ( |2.aD one readily finds 



"(X, t) = U N (x) TJ i U b ■ U *N ( X b)} ReLators.TV (*) + Yl /<Sernal,AT RxLnal^ (*) ) ( 70 ) 

N 1 6=1 « ) 

where the label "external" explicitly refers to agents acting upon the system from outside. Two kinds 
of such external actions are considered in this article: those due to GWs and those due to a calibration 
hammer stroke signal. Specifying x = x a in the lhs of ( |70| ) and multiply on either side by n a , the 
following is readily found: 



Ua(t) = < 



external 



(*) + Z> f K ab (t - t') [ z b (t') - u b (t')] dt' (71.a) 

6=1 J ° 



%a(t) = C temal (t) " tf a [z a (t) - u a (t)\ , a = 1, . . . , J (71.b) 
where % = M b /M , n ^ xternal (t) = n a ■ u cxternal (x a , t) , 

u e x ternal (X) t) = £ £ ^ Ujy(x) g^^t) ( 7 2) 

a N 

and 

K ab (t) = Of, ^n 1 [ n 6' u Af( x 6)] [n a -UAr(x a )] sin u N t (73) 

N 

The following bare sphere responses to GWs and hammer strokes (equations! @j and (^), respec- 
tively) can be calculated by direct substitution. The results as best presented as Laplace transform 
domain functions: 

u™(s) = £ (±^^)Y lm (n a )9 (lm) (s), a = l,...,J (74a) 

Z=0 and 2 \n=l nl / 

m=— 

.stroke^ = _^_A_|^ (i?) |2^ (na . no); a = 1, . . . , J (74.b) 

nl S +UJ nl 
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where Y\ m are spherical harmonics and Pi Legendre polynomials [15|. The calculation of the Laplace 
transform of the kernel matrix ([7^) is likewise immediate: 



K ab (s) = £ 



N 



s 2 +uj 2 n 



[n b -u* N (x b )] [n a -U7v(x a )] 



(75) 



Given that (see 24] for full details) 



u„; m (x) = A ni (r) Y lm (0, <p) n- B nl (r) in xLY lm (8,(p) (76) 

and that the spheroidal frequencies uj n i are 2Z+l-fold degenerate, ( [75| ) can be easily summed over the 
degeneracy index m, to obtain 



nl 



J2 Y^b)Y lm (n a ) 

m=—l 



or, equivalently, 



(77) 



A- \A nl (R)f «±1 



nl 



S z + LU. 



nl 



47T 



Pi(n a -n b ) 



where use has been made of the summation formula for the spherical harmonics [15] 



(78) 



21 + 1 

E Y l*m( n b) ^m(n a ) = — — ^( n a-n fe ) 

m=— I 



Air 



and where Pi is a Legendre polynomial: 



(79) 



(80) 



A. 2 System response algebra 

From equation (^3|), i.e., 



E 

6=1 



2„2 



ft 2 S 



(nJ) 



(s 2 + 02)(s2+ ^ ) X o6 



i GW (*) 

s 2 + ft 2 



s 2 + n 2 



(81) 

we must first isolate cjb(s), then find inverse Laplace transforms to revert to time domain quantities. 
Substituting the values of u% w (s) and ^ w (s) from ( frTal) and (||) into (||) we find 



q a (s)= E ^ m) (^)5 (/m) (^) 



a = 1 J 



(82) 



Z = and 2 
m = — Z, I 



where 



s 2 + ft 2 



e 2 + 



n=l 



E 

f>=l 



^afc + V E 



ft 2 S 



2„2 



(ni) 



nl 



(s 2 + n 2 )(s 2 + oj nl )< 



(83) 

Now, using the convolution theorem of Laplace transforms, we see that the time domain version of 
equation (82) is 
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9a(*)=E f ^ m) (t-t')g {lm) (t')dt' 



,J 



(84) 



where $>a m \t) is the inverse Laplace transform of (|33"1), The inverse Laplace transform of <& a lm \s) can 
be expediently calculated by the residue theorem through the formula 



<fr a lm \t) = 2iri E {residues of <l^ m ^(s) e st at its poles in complex s — planej 



(85) 



Clearly thus, the poles of &a m \s) must be determined in the first place. It is immediately clear 
from equation (^) that there are no poles at either s = 0, or s = ±zO, or s = ±iui n i, for there are exactly 
compensated infinities at these locations. The only possible poles lie at those values of s for which the 
matrix in square brackets in ( p3| ) is not invertible, and these of course correspond to the zeroes of its 
determinant, i.e., 



A(s) = det 



b~ab + V ■ 



K ab (s) 







poles 



(86) 



There are infinitely many roots for equation (|8q), but analytic expressions cannot be found for 
them. Perturbative approximations in terms of the small parameter rj will thus be applied instead. It 
is assumed that 



^ = w„ i (87) 
for a fixed multipole harmonic {uqIq}. Equation (|86|) can then be recast in the more convenient form 



A(s) = det 



0~ab + V 



(s 2 + n 2 ) 



2 x a ? l0) + v E 



n 2 s 2 



(nl) 



{s 2 + n 2 )(s 2 + co 2 nl ) Xab 







Since n is a small parameter, the denominators of the fractions in the different terms in square 
brackets in ( p8| ) must be quantities of order r\ at the root locations for the determinant to vanish at 
them. A distinction however arises depending on whether s 2 is close to — VL 2 or to the other —co^i- 
There are accordingly two categories of roots, more precisely: 



4 



,2 
s nl 



-n 2 (i + xi r/ 1/2 + xi v + ■ ■ •) (n = ^Wo) 

-u 2 nl (l + &f"> r? + r, 2 + . . .) (nl + n lo) 



(89.b) 



The coefficients yi, xi,... and 6i , &j can be calculated recursively, starting form the first, 

2 

by substitution of the corresponding series expansions into equation ( pq ) . The lowest order terms are 
easily seen to be given by 



det 



c (no Jo) 

d «6 - -2" X„ 6 

2 



(90) 



and 



det 



9 2 



^nl h (nl) , (nl) 



nl 



(91) 



respectively. Both equations ([H]) and ( |9l| ) are algebraic eigenvalue equations. As shown in ap- 



pendix A.3 , the matrix xlfe has at most (21 + 1) non-null positive eigenvalues — all the rest up to 
J are identically zero. 
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As a final step we must evaluate (|8q) . This is accomplished by standard textbook techniques (see 
e.g. the algebra is quite straightforward but rather lengthy, and I shall not delve into its details 

here, but quote only the most interesting results. It appears that the dominant contribution to <&a \t) 



comes from the poles at the locations (89. a), whereas all other poles only contribute as higher order 
corrections; generically, $>a m \t) is seen to have the form 



where 



^ m Ht) oc £ ( si nu; c+ t - sinc^t) «5„ + O(0) 

Cc^o 



(92) 



(>' \ l±J 2 l±± \A n0l0 (R)\CaV 1/2 \ +O( V ) 



a = 1 J 



(93) 



In Laplace domain one has, 



1/2 E ( 



s 2 + uj 2 c+ 



r 1 -(*+<*-) 



(94) 



Detailed calculation of the residues pq ] yield equation (Pq), which must be evaluated for each 
particular tuning and resonator distribution, as described in section 0. 



A. 3 Eigenvalue properties 

This Appendix presents a few important properties of the matrix P/(n a -n;,) for arbitrary I and resonator 
locations n a (a=l,. . . ,J) which are useful for detailed system resonance characterisation. 
Recall the summation formula for spherical harmonics [15]: 



21 + 1 

E Y lm( n a)Ylm(nb) = A _ Pl^a'W) , 

m=—l 



47T 



a,b = 1, . . . , J 



where Pi is a Legendre polynomial 



2 x y 



(95) 



(96) 



To ease the notation I shall use the symbol V\ to mean the entire JxJ matrix P;(n a -n(,), and 
introduce Dirac kets \m) for the column J- vectors 



771) = 



/ 4vr 
2Z + 1 



^m(ni) 



777 = —I, . . . , I 



(97) 



These kets are not normalised; in terms of them equation ( p5| ) can be rewritten in the more compact 
form 



Pi = ^2 \m)(m\ 



(98) 



m=—l 



Equation ( |9q ) indicates that the rank of the matrix Vi cannot exceed (21 + 1), as there are only 
(21 + 1) kets \m). So, if J > (21 + 1) then it has at least (J — 21 — 1) identically null eigenvalues — there 
can be more if some of the n a 's are parallel, as this causes rows (or columns) of Vi to be repeated. 

We now prove that the non-null eigenvalues of V\ are positive. Clearly, a regular eigenvector, \<p), 
say, of Vi will be a linear combination of the kets |m): 
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I 

vi 10) = c 2 10) , 10) = E 0™ H (") 

m=— I 

where ( 2 is the corresponding eigenvalue, having a positive value, as we now prove. If the second ( p9[ ) 
is substituted into the first then it is immediately seen that 

i 

E (C 2 <W - (m|m')) 0m' = (100) 
m'=—l 

which admits non-trivial solutions if and only if 

det (C 2 5 mm > - (m|m')) = ( 101 ) 

In other words, ( 2 are the eigenvalues of the (21 + l)x(2Z + 1) matrix (m\m'), which is positive 
definite because so is the "scalar product" (0|0'). All of them are therefore strictly positive. 
Finally, since the trace is an invariant property of a matrix, and 

J J 
trace(^) = E ^K-n a ) = E 1 = J ( 102 ) 

a=l a=l 

we see that the eigenvalues Q 2 add up to J: 

J 

trace(^) = ^C 2 ^( a 2 = J (103) 

a=l C a ^0 
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Figure 2: The three distinct eigenvalues £ m (m = 0,1,2) as functions of the distance of the resonator 
parallel's co-latitude a relative to the axis of symmetry of the distribution, cf. equation ( |39.a| ). 
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Figure 3: To the left, the pentagonal hexacontahedron shape. Certain faces are marked to indicate 
resonator positions in a specific proposal — see text — as follows: a square for resonators tuned to the 
first quadrupole frequency, a triangle for the second, and a star for the monopole. On the right we see 
the (pentagonal) face of the polyhedron. A few details about it: the confluence point of the dotted lines 
at the centre is the tangency point of the inscribed sphere to the PHC; the labeled angles have values 
a = 61.863°, [3 = 87.205°; the angles at the T-vertices are all equal, and their value is 118.1366°, while 
the angle at P is 67.4536°; the ratio of a long edge (e.g. PT\) to a short one (e.g. T1T2) is 1.74985, 
and the radius of the inscribed sphere is twice the long edge of the pentagon, R = 2PT±. 
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Figure 4: Compared line spectrum of a coupled TIGA and a PHC resonator layout in an ideally 
spherical system. The weakly coupled central frequency in the TIGA is drawn dashed. The frequency 
pair is 5-fold degenerate for this layout, while the two outer pairs of the PHC are doubly degenerate 
each, and the inner pair is non-degenerate. Units in abscissas are r] l / 2 Vl, and the central value, labeled 
0.0, corresponds to Q. 
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Figure 5: Simulated response of a PHC to a hammer stroke: the time series and their respective spectra, 
both for direct resonator readouts and mode channels. Note that while the former are not simple beats, 
the latter are. 
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Figure 6: The frequency spectrum of the TIG A distribution as resonators are progressively added from 
none to 6. Continuous lines correspond to measured values, and dashed lines correspond to their r/ 1 / 2 
theoretical estimates with equation (56). 
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